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Abstract. We are concerned with the inviscid limit of the Navier-Stokes equa- 
tions to the Euler equations in M 3 . We first observe that a pathwise Kolmogorov 
hypothesis implies the uniform boundedness of the a t,l -order fractional derivative 
of the velocity for some a > in the space variables in L 2 , which is independent 
of the viscosity /i > 0. Then it is shown that this key observation yields the 
L 2 -equicontinuity in the time and the uniform bound in L q , for some q > 2, of 
the velocity independent of \i > 0. These results lead to the strong convergence 
of solutions of the Navier-Stokes equations to a solution of the Euler equations in 
R 3 . We also consider passive scalars coupled to the incompressible Navier-Stokes 
equations and, in this case, find the weak-star convergence for the passive scalars 
with a limit in the form of a Young measure (pdf depending on space and time). 
Not only do we offer a framework for mathematical existence theories, but also 
we offer a framework for the interpretation of numerical solutions through the 
identification of a function space in which convergence should take place, with the 
bounds that are independent of fi > 0, that is in the high Reynolds number limit. 



I. Introduction 

Consider the incompressible Navier-Stokes equations in M?\ 

d t u + V ■ (u <g> u) + Vp = fiAu + f , 
V-u = 0, 



with Cauchy data: 

u| t=0 = u (x), (1.2) 

where u is the fluid velocity, p is the pressure, /i > is the viscosity, V is the gradient 
with respect to the space variable x 6 M 3 , u <g> u = (uiUj) is the 3x3 matrix for 
u = (ui, u 2 , u 3 ), and f = f (t, x) is a given external force. 

The global existence theory for the Cauchy problem f ll.ip - fll.2p was first estab- 
lished by J. Leray [HCL5JCL6]; also see Hopf [IT], Temam [27], J.-L. Lions [19], P.-L. 



Date: November 2, 2011. 

2010 Mathematics Subject Classification. 76D05, 35Q30, 76D09, 76F02, 65M12. 

Key words and phrases. Inviscid limit, convergence, Navier-Stokes equations, vanishing viscos- 
ity, Kolmogorov's hypotheses, turbulence, equicontinuity, existence, weak solutions, Euler equa- 
tions, numerical convergence. 

1 



2 



GUI-QIANG CHEN AND JAMES GLIMM 



Lions [20], and the references cited therein. For clarity of presentation, we focus on 
periodic solutions with period T P = [-P/2, P/2] 3 cR 3 ,P> 0, that is, 

u^x + Pe,;) = u"(t,x) 

with (ej)i<j<3 the canonic basis in M 3 . Other cases can be analyzed correspondingly. 
We always assume that f G Lf oc (0, oo; L 2 (T P )) is periodic in x with period Tp. 

Theorem 1.1. Let Uq G L 2 (Tp) be periodic in x G IR 3 with period Tp. Then, 
for any T > 0, there exists a periodic weak solution u M = u M (t,x) with period Tp, 
along with a corresponding periodic pressure field p M (£, x) , of (jl.ip - (jl.2p such that 
the equations in (II. ip hold in the sense of distributions in IR^ := [0,T) x IR 3 , and 
the following properties hold: 

u» G L 2 (0, T; H 1 ) H C([0, T]; L 2 ) PI C([0, T]; L* 1 ), (1.3) 

t u" G L^O^PT 1 ) + (L S2 (0,T;^ -1 '^) nL 9 (0,T;L r )), (1.4) 

p M G L 2 ((0,T) x T P ) + L S2 (0,T;L^), (1.5) 
Vp M G L 2 (fd,T-H- 1 ) + P ? (0,T;Z7), (1.6) 

where 1 < si < 2, 1 < s 2 < oo, 1 < q < 2, and r = 2 (2q-i) >' an< ^' ^ n addition, 

<9 t (i|uf ) + V- (u^(i|u^| 2 +^)) +/i|Vu^| 2 -//A(^£) < f -u^ (1.7) 

m t/ie sense o/ distributions in IR^. 

In Theorem 1.1, v G C([0, T]; L 2 (T P )) means that v G L°°(0, T; L 2 (T P )) and v 
is continuous in t with values in L 2 (T P ) endowed with the weak topology. Some 
further a priori estimates and properties of solutions to the Navier-Stokes equations 
fll.lD can be found in El [20] and the references cited therein. 

In contrast, less is known regarding the existence theory for the Euler equations 
(defined with \i = in (jl.lft ). For the compressible case, the analysis of [TUJ |22] 
gives convergence subsequences, but to a very weak limit as a measure- valued vector 
function. Moreover, this limit is not shown to satisfy the original equations, in 
that the interchange of limits with nonlinear terms in the equations is not justified 
in this analysis. On this basis, we state that the existence of solutions for the 
Euler equations in M 3 is open as is the convergence of the inviscid limit from the 
Navier-Stokes to the Euler equations. On the other hand, the Euler equations are 
fundamental for turbulence; see Constantin [1] and the references cited therein. 

The purpose of this paper is to establish a framework for the existence theory for 
the Euler equations. We introduce a physically well accepted pathwise hypothesis by 
Kolmogorov [I2l ll3|. Assumption (KJ^l). One of our main contributions in this paper 
is our key observation that Assumption (K41), even a weaker version Assumption 
(K41w), assures sufficient regularity of the solutions to the Navier-Stokes equations 
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with an external force that convergence through a subsequence to solutions of the 
Euler equations is guaranteed. 

As is well-known, there are two types of turbulence: driven turbulence by a forcing 
function and transient turbulence by (strong) initial conditions. The same flow (in a 
turbulent wind tunnel or turbulent flow in a pipe) could be of either type depending 
on how the system is modeled: If the pipe is considered in isolation, the flow might 
be forced; if the force is a flow connection from a reservoir to the wind tunnel or 
pipe, and the reservoir is part of the model, then the turbulence arises from initial 
conditions, and the turbulence will die out when the reservoir is exhausted and no 
longer drives the flow. Thus, the distinction between the two (transient and driven 
turbulence) is to some extent a matter of points of view and of modeling convenience. 
In this paper, our framework is focused on transient turbulence, though the forced 
turbulence is also included. For some recent developments in the mathematical 
study of energy dissipation in body-forced turbulence, see Constantin-Doering [5], 
Doering-Foias [8], and the references cited therein. 

Whenever a mathematically precise formulation of a physically precise idea is 
attempted, ambiguities may arise. In the present case, the assumptions concern 
the rate e of dissipation of kinetic energy, with a fundamental definition which 
depends on the viscosity, as the ultimate source of energy dissipation. The essence of 
Assumption (K41), as stated in a physical language, is that in the inertial range, the 
energy is transferred from larger to smaller length scales (eddies) in a manner which 
is independent of viscosity, because all aspects of the inertial range are assumed to 
be independent of viscosity. The energy dissipation occurs only at smaller scales, 
i.e., below the Kolmogorov scale, and even there, it is limited by the energy arriving 
at these scales through the energy cascade. Given that the energy dissipation rate is 
supposed to be independent of viscosity, the question remains as to which statistical 
ensemble or in what norm or function space to express this property. Here the 
physical literature is not a good guide, and we introduce the assumptions sufficient 
to allow our proofs to go forward. The distinction among time averages, spatial 
averages at a fixed time, space-time averages and ensemble averages relates to the 
well known ergodic hypothesis and is out of the scope of the present paper. 

Mathematically, our result has the status of an informed conjecture and math- 
ematically rigorous consequences of this conjecture. Numerical analysts may find 
the framework useful, in view of the many difficulties involved in assessing conver- 
gence of numerical simulations of turbulent and turbulent mixing flows. We expect 
that many physicists will probably accept the conclusions as being correct, even if 
unproven mathematically. There has been some discussion regarding the exponent 
5/3 which occurs in Assumption (K^l). We note that the main results (if not the 
detailed estimates) are not sensitive to this specific number, and corrections (as con- 
ventionally understood) to it due to intermittency should not affect our result. In 
fact, our rigorous argument works for an even weaker version, Assumption (KJ^lw), 
for any (3 > 0. 
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Not only do we offer a framework for mathematical existence theories, but also 
we offer a framework for the interpretation of numerical solutions of (II. ip . Only 
for very modest problems and with the largest computers can converged solutions 
of (II -ip be achieved. These solutions are called direct numerical solutions (DNS). 
For most solutions of interest to science or engineering, the large eddy simulations 
(LES) or Reynolds Averaged Navier-Stokes (RANS) simulations are required. We 
discuss here the more accurate LES methodology. Briefly, it employs a numerical 
grid which will resolve some but not all of the turbulent eddies. The smallest of 
those, below the level of the grid spacing, are not resolved. However, either (a) the 
equations in (II. ip are modified with additional subgrid scale (SGS) terms to model 
the influence of the unresolved scales on those that are being computed or (b) the 
numerical algorithm is modified in some manner to accomplish this effect in some 
other way. The present article contributes to this analysis through the introduction 
of a function space in which convergence should take place, with bounds that are 
independent of /i, that is in the high Reynolds number limit. 

Because of the common occurrence of high Reynolds numbers in flows of practical 
and scientific interest and the need to perform LES simulations to achieve scientific 
understanding and engineering designs, we observe that existence theories for the 
Euler equations are relevant to the mathematical theories of numerical analysis. 



2. The Kolmogorov Hypotheses (1941) 



By the definition of weak solutions, we know that, for any T > 0, there exists 
Ct > independent of /i such that 

\\ U>1 ~ " Q lll=°(0,T;L 2 (Tp)) + II V/^ U '1l-L 2 ([0,T)xTp) 
< C T (||u - u||| 2 ( Tp ) + ||f - f ||z,2([ 0) T)xTp))> 



where 



1 



u(t) = — u ^,x)dx=-— / u (x)dx+ / f(s)ds 

l^Pl Jfp \ & P\ JTp JO 

and 

f(f) = J-, / f(f,x)dx 

are independent of /i. Without loss of generality, we assume that the mean velocity 
and pressure are zero, and interpret as the fluctuating velocity. Then the total 
energy S(t) per unit mass at time t for isotropic turbulence is: 

Here E(t, k), k = |k|, is the energy wavenumber spectrum, q(t, k) can be interpreted 
as the density of contributions in wavenumber space to the total energy, which is 
sometimes called the spectral density, and k = (ki,k 2 ,k 3 ) = ^(^i,n 2 ,n 3 ) e M 3 , 
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with rij = 0, ±1, ±2, • • • , and j = 1, 2, 3, is the discrete wavevector in the Fourier 
transform: 

u(t,k) = -^ / U (t,x)e- ik ' x rfx 

of the velocity u(t, x) in the x- variable. Then 

u(t,x) = ^u(t,k)e lk ' x 

k 

Kolmogorov's two assumptions in his description of isotropic turbulence in Kol- 
mogorov [T21 [32] (also see McComb [21]) are: 

(i) At sufficiently high wavenumbers, the energy spectrum E(t,k), can depend 
only on the fluid viscosity /z, the dissipation rate e, and the wavenumber k 
itself. 

(ii) E(t,k) should become independent of the viscosity as the Reynolds number 
tends to infinity: 

E(t,k)^ae 2/3 k- 5 / 3 (2.3) 

in the limit of infinite Reynolds number, where a may depend on t, but is 
independent of k and e. 

For general turbulence, the energy wavenumber spectrum E(t, k) in i \2.3\\ may be 
replaced by E(t, k) = E(t,k,(p,8) in the spherical coordinates (k,(p,8), < < 7r, 
< 6* < 27r, in the k-space, but it should be in the same asymptotics as in (I2.3P for 
sufficiently high wavenumber k = |k|. 

For the remainder of this paper, we assume these two hypotheses, which we in- 
terpret in mathematical terms as a pathwise Kolmogorov hypothesis: 

Assumption (K41): For any T > 0, there exists Ct > and (sufficiently 
large) depending on uq and f but independent of the viscosity [i such that, for k = 
|k| ^ k.%, 

T 

E(t,k)dt<C T k- 5/3 . (2.4) 

We should point out that the constant Ct = Cj^uq, f) depends on the maximum 
time T, the initial data u , and the external force f, but is independent of the 
viscosity fi > 0. For given the maximum time T, the dependence on Uo and f 
of the constant Ct is from the control of the upper bound of J Q T e 2 ^(s)ds, with 
e = e(t) being the rate of energy dissipation, by the initial data and the external 
force, more specifically by ||uo ||i 2 (x^) an d ||f ||l 2 ([o,t)xt p )- The dimensionless version 
of the bound follows from (12. 3p . 



For our analysis, the following weaker version of Assumption (K41) is sufficient: 
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Assumption (K41w): For any T > 0, there exists Ct = Ct > and 
(sufficiently large) depending on u and f but independent of the viscosity jj, such 
that, for k = |k| > k*, 



A mathematical proof of Assumption (K41) may well depend on developing a 
mathematical version of the renormalization group. The renormalization group has 
proved to be very powerful in theoretical physics calculations. The basic idea is to 
integrate differentially a segment of the problem, say from k to k — d\t, following 
by a rescaling of all variables, so that all variables are redefined to be integrated 
through k only. The integration is quite a messy operation, but in the rescaling, 
it is important to examine the scaling dimensions of all terms. Most of them get 
smaller under rescaling and are called inessential. The few that are preserve their 
magnitude are called essential and serve to define the key parameters of the problem. 
These must be reset (renormalized) to their observed (interaction) values after the 
integration. Since the method is generally applied to self-similar problems, the result 
of integration does not change the form of the equations. That is, it looks the same 
at k as at k — d\t rescaled back to k. Thus, one has a kind of group operation, if the 
steps are discrete, and a differential equation, if the steps, as the present notation 
suggests, are infinitesimal. In either case, the desired solution is a fixed point. 



3. L 2 -Equicontinuity of the Velocity in the Space Variables, 
Independent of the Viscosity 

In this section we show that the pathwise Kolmogorov hypothesis, Assumption 
(KJ^lw), implies a uniform bound of the velocity u^(t, x) in L 2 (0, T; H a (Tp)) for 
any a G (0,/3/2), especially the uniform equicontinuity of the velocity u M (t, x) in 
the space variables in L 2 ([0,T) x Tp), independent of fj, > 0. 

Proposition 3.1. Under Assumption (K41w) ; for any T e (0,oo) 7 there exists 
Mt > depending on u , f , and T > 0, but independent of ji > ; such that 




for some (3 > 0. 



(2.5) 



U M ||L2( 0jT;//a ( Tp )) < M T < 00, 



(3.1) 



where a E (0,(3/2). 
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Proof. Using the definition of fractional derivatives via the Fourier transform, the 
Parseval identity, and Assumption (K41w), i.e., f)2.5p . we have 

/ \D°rt l (t,x)\ 2 dxdt 

'o JJp 

< 



G x [ (^|k| 2o |u>(t,k)| 2 )dt 
Jo k 

^ f { E Ikl^lu^kJI^ + d f ( Ik| 2 »(t,k)| 2 ) 



0<|k|<fc» " u |k|>fc» 
r-T 

\2a-3-f3 



< c ^ k * a I ( E iu>(t,k)i 2 )^+c 2 ^ |k 

0<|k|<fc* |k|>fe» 

< ^* 2Q f (Ei^^k)! 2 )^ +^E 

^° k fc>fc* 

< C^ 20 / f |u"(t, x)| 2 dxdt + Cik 2 *' 13 

Jo Jj P 

<c 5 k 2a r( [ | Uo (x)| 2 rfx+ [ T [ \f(t,x)\ 2 dxdt)+c,kl a -p 

V JT P JO Jf P ' 

< M(T, K,a) < oo, 

since a < (3/2, which Cj,j = 1, • • ■ ,5, are the constants independent of /i. This 
completes the proof. □ 

Proposition 3.1 directly yields the uniform equicontinuity of u M (t,x) in x in 
L 2 ([0,T) x T P ) independent of /i > 0. 

4. L 2 -Equicontinuity of the Velocity in the Time- Variable, 
Independent of the Viscosity 

In this section, we show that Proposition 3.1 implies the uniform equicontinuity 
of the velocity in the time variable t > in L 2 , independent of \i > 0. 

Proposition 4.1. For any T > 0, there exists > depending on u ,f, and 
T > 0, but independent of > 0, such that, for all small At > 0, 

rT-At r 

|u M (t + At,x) - u M (t,x)| 2 dxdt < M T (At)5T^ as At -> 0. 



JTp 

(4.1) 

Proof. For simplicity, we drop the superscript fi > of u M in the proof. 
Fix At > 0. For te[0,T- At], set 

w(t,-) = u(t + At,-) -u(t,-). 
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Then, for any ip(t,x) G C°°([0,T) x Tp) that is periodic in x 6 R 3 with period Tp, 
we have 



w(i, x) • (p(t, x) (ix 



rt+At r- 

i / <9 s u(s,x) • x) dxds 

Jt JTp 
t+At p 

/ (u ® u)(s, x) : V x) dxds 
+ / / p(s, x) V • cp(t, x) dxds 

/■t+At p 

—fj, / / Vu(s, x) : V<p(t, x) (beds 
Jt Jip 



rt+At r 

+ / / f(s,x) • (p(t,x)dxds, (4.2) 



where = (d Xj <pi) is the 3x3 matrix, and A : 5 is the matrix product ^ . dy-ify 
for A = (a,ij) and 5 = 

By approximation, equality ( 14. 2 p still holds for 



cp G L oo (0,T;H 1 (T P )) H C([0, T]; L 2 (T P )) nL°°([0,T) x T P ). 



Choose 



<P = V 5 (*, x ) : = (j<5* w)(t,x) = / j 5 (y)w(t,x-y)c?y G 



for 6 > 0, which is periodic in x G M 3 with period Tp, where .7,5 (x) = ) > is 
a standard mollifier with j G C^°(1R 3 ) and J j'(x)c?x = 1. Then 



V V(t,x) = 0, (4.3) 



since V • w(i, x) = 0. 
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Integrating (14.21) in t over [0, T — At) with cp = tp s (t, x) and using (14.31) . we have 



T—At p 

/ |w(t, x)\ 2 dxdt 

Jf P 

T—At rt+At 



pi — 1st pt-\-i\t p 

/ / / (u®u)(s,x) : Vip s (t,x)dxdsdt 

J0 Jt JT P 

pT-At r pt+At 

-fi / / Vu(s,x) : Vy?*(t,x) dsdxdt 

Jo Jj p Jt 

pT-At p 

+ / / w(t, x) • (w(i, x) - (p s (t, x)) dxdt 
Jo Jl P 



+ 



T—At r t+At 



t JT 



f(s,x) • (p s (t,x)dx.dsdt 



4 + 4 + 4 + 4- 



(4.4) 



Notice that, for any x e Tp, 



|V</(t,x)| < 



1 



5 



< 



< 



f s (x-y)w(t,y)dy 

x-y|<<5 

§f 2 { J |j"(y)l 2 rfy) 1/2 ||w(t,-)llL 2( T P ) 
c 



§5/2- 



Here and hereafter, we use C > as a universal constant independent of /i > 0. 
Then 



s CTAt 2 CTAt 2 CAt 

l J ll - £5/2 II U (^-)IIl2(T p ) < J5/2 II U 0|Il2 (Tp) < 



(4.5) 
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Furthermore, we have 

pT—At r r t+At 



< 



JTp Jt 

T—At I- , r t+At 



pi — lAt p pt + lAt 

\J$\ < / / fi\Vu(s,x.)\\V(p s (t,x)\dsdxdt 

Jo Jjp Jt 

v //I|Vu(s,x)|dsj v //I|V(p <5 (t,x)| dxdt 

pl-l\l p pt+t\t 

< ( / / ( / V^|Vu(s,x)|ds) dxdt 

JO J Tp J t 



T—At 



/j,\Vcp d (t,x.)\ 2 dxdt 

T—At r rt+At 



1/2 



pi —LAt p pt+lAt 

< (At) 1/2 ( / / / /i|Vu(s,x)| 2 dsdxcft 

^ Jo Jt p Jt 

x( f f n( f j 5 (y)|Vw(t,x-y)|dy) 2 rfxdt) 1/2 
v Jo JTp J\y\<S J 

< CAtL f \j 5 (y)\ 2 dy [ [ [ fi\ Vw(t, x - y)\ 2 dxdtdy 

V Jm. 3 J\y\<&Jo Jj P 



1/2 



r |y|<<5. 

< CAt, (4.6) 
where we have used 

Hv7^Vu|lia([0,T)>cTp) < m t 
from (11. 7p for the weak solutions in Theorem 1.1 with > independent of \x. 
On the other hand, for 

pT—At p 

J 3 5 := / / w(t,x)(w(t,x) - (p s (t,x))dxdt, 

Jo Jjp 



we find 



| J|| = / / w(t,x)( / j 5 (x-y)(w(t,x)-w(t,y))dyWt 

J0 JTp X JR 3 J 

pT—At p p 

< / / j(y)|w(t,x)||w(t,x)-w(t,x-5y)|dydxdt 

•/() Jfp J\y\<l 

< / 3(y)( I I |w(t,x)||w(t,x)-w(t,x-5y)|dxrft)rfy^ 

J\y\<l y Jo Jjp / 

< C [ j(y)( f [ |w(t,x)-w(t,x-5y)| 2 rfxrft) 1/2 rfy^ 

J\y\<l y Jo Jt p ' 

< C5 a [ j(y)\y\ a dy 

jr 3 

< C 3 5 a . (4.7) 
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Here we have used the fact that ||w(£, -)I|l 2 (Tp) < C and 

|w(t, x) — w(£, x — 5y)\ 2 dx.dt 



o J J 



< 



/ T (^|w(t,k)| 2 (l-e^) 2 )^ 
k 

C5 2a \y\ 2a f T (^|k| 2a |w(t,k)| 2 )dt 
k 

<C5 2a \y\ 2a [ [ \D^(t,x)\ 2 dxdt 

Jo Jf P 

< C6 2a \y\ 2a . 



Moreover, we have 



|jf| = CAt\\f ||L 2 ([0,T)xTp)||V 5 ' 5 |U 2 ((0,r-At)xTp) 

< CAt||w|| L 2(( 0jT _ At ) xTp) 

< CAt||u|| L 2(( 0jT ) xTp ) 

< C A M. (4.8) 
Combining fr4T4]) -f H~Tj) with (Q|) . we have 



i-T-M p 
JO Jj f 



,T-At r a. 

|w(t,x)| 2 dxdt < inf {C x + (C 2 + C 4 ) At + C 3 5 a } 



5>0 L °5 5 / 

Choose 

5C; 



We conclude 

eT-At 



j I \w(t,x.)\ 2 dxdt<C(At)^. (4.9) 

This completes the proof. □ 

As a direct corollary, we have 

Proposition 4.2. Under Assumption (K41w) ; for any T G (0, oo), there exists 
Mt > depending on uq, f , and T > 0, but independent of fi > 0, such that 



[ [ llWftx^dtdx < M T < oo, (4.10) 
Jo Jf P 



for some a > 0. 
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Combining Proposition 3.1 with Proposition 4.2, we have 

||u m ||h-([o,t)xTp) < M T < oo (4.11) 

for some a > 0. Then, by the Sobolev imbedding theorem, we have 

Proposition 4.3. Under Assumption (K41w) ; for any T G (0,oo) ; there exists 
Mt > depending on uq, f , and T > 0, but independent of \x > 0, such that 

||u M || L 2nL9([o,T)xTp) < M T < oo (4.12) 

for some q > 2. 



5. Inviscid Limit 

In this section we show the inviscid limit from the Navier-Stokes to the Euler 
equations under Assumption (KJ^lw). 

Theorem 5.1. The pathwise Kolmogorov hypothesis, Assumption (K41w) ; implies 
the strong compactness in L 2 fl L q ([0,T) x Tp), for some q > 2, of the solutions 
u M (t, x) of the Navier- Stoke equations in K 3 when the viscosity [i tends to zero. 
That is, there exist a subsequence (still denoted) u^(t,x) and a function u e L 2 fl 
L q ([0,T) x Tp) with a corresponding pressure function p such that 

u M (t, x) — > u(t, x) a.e. as \i — > 0, 

and u(t, x) is a weak solution of the incompressible Euler equations with Cauchy 
data uo(x) along with the corresponding pressure p. Furthermore, 

[ |u(t,x)| 2 dx< f |u (x)| 2 dx + 2 [ [ u(s,x) ■ f(s,x)rfxrfs. (5.1) 

JTp JTp JO JTp 



Proof. Propositions 3.1 and 4.1 imply the L 2 -equicontinuity of u M (t, x) in (t, x) G 
[0, T) x Tp, independent of /x > 0. This yields that there exists a subsequence (still 
denoted) u M (t, x) and a function u(t, x) G L 2 such that 

u^it, x) u(t, x) in L 2 as /i -> 0, 

which implies that 

u M (t, x) — > u(t,x) a.e. as — > 0, 
(u^ <g> u M )(t,x) (u (g) u)(t,x) a.e. as /i 0, (5.2) 

and 

V-u(t,x)=0 (5.3) 
in the sense of distributions in 1R T . 
From Proposition I4.3[ we have 

u G L 2 nL 9 ([0,T) x Tp), 
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and 

u M (£, x) ->■ u(f, x) in L 2 n L 9 ([0, T) xT P ) as/i^O 
for some q > 2. 

Furthermore, for any e C^°(Ry) with V • (f = 0, we multiply <£> both sides (11. ip 
and integrate over Rj. to obtain 

[ [ (u>V t + (u^ ® u**) : Vy> + • f <p) dxdt + / u o (x)^(0, x) dx 

= -n I I u M ■ A<p dxdt. 
Jo Jr 3 

Letting \i — > 0, we conclude that 

/ f (u<p t + (u (8) u) : V<p + (u • f) <p) tZxdi + / u (x)y>(0, x) dx = 0. (5.4) 

</0 ./R 3 </R 3 

Combining (15 .4p with ( 15. 2 p yields that u(t, x) is a weak solution of the Euler equa- 
tions (II. ip (/i = 0) with Cauchy data (ll.2p along with the corresponding pressure 
p. 

Integrating (11.7]) over [0,T) x Tp yields 



ft 

|u^(t,x)| 2 dx< / |u (x)| 2 dx + 2 / / u%s,x) ■ f(a,x)dxds. 

Tp •/• ./() JTp 

Letting \i — > 0, we have 

/ |u(t,x)| 2 dx< / |u (x)| 2 dx + 2 / / u(s,x) ■ f(s,x)t2xds. 

Jt p JTp JO JTp 

This completes the proof. □ 

Remark 5.1. The question whether the energy inequality in ( 15. ip becomes identical 
is exactly Onsager's conjecture [23]. It states that weak solutions of the Euler 
equations for incompressible fluids in R 3 conserve energy only if they have a certain 
minimal smoothness of the order of 1/3 fractional derivatives in x e R 3 and that they 
dissipate energy if they are rougher. Indeed, in Cheskidov-Constantin-Friedlander- 
Shvydkoy [2J, it is proved that energy is conserved when u(t, x) is in the Besov space 

1/3 

of tempered distributions -B 2 c (N)- This is a space in which the Holder exponent 
is exactly 1/3, but the slightly better regularity is encoded in the summability 
condition. We show in Proposition 3.1 that Assumption (KJ^lw) implies the uniform 
boundedness in the norm of L 2 (0, T; H a (Tp)) for some a > for the solutions to 
the Navier-Stokes equations, independent of the viscosity fi > 0. From the point 
of view of physics, strengthening the Kolmogorov exponent 5/3 (i.e. /3 = 2/3), as 
would be needed to obtain a stronger a = 1/3 bound, may result from some theory 
of intermittency. 



14 gui-qiang chen and james glimm 

6. Incompressible Navier-Stokes Equations with Passive Scalar 

Fields 

Consider the incompressible Navier-Stokes equations with passive scalar fields in 
IR 3 : 

d t u + V • (u <g> u) + Vp = yU Au + f , 

dtXi + V • (x<u) = /x»Axi, « = 1, • • • (6.1) 
V-u = 0, 
with Cauchy functions: 

(u,Xi,--- ,Xi)\t=o = (u ,Xio, ••• ,X/o)(x) (6.2) 

that are periodic in x e M 3 with period T P , where x«o( x ) > with XliXio( x ) = 1> 
and f G £ 2 OC (0, oo; L 2 (Tp)) is a given external force. 

These equations are interpreted as describing concentrations Xii i < I °f minority 
species, such as dilute chemical reagents in a majority carrier species xii m an 
approximation that the minority species do not influence the bulk fluid density nor 
the bulk fluid viscosity. Thus, they do not interact with (influence) the bulk fluid, 
but are passively transported by it; hence the name passive scalars. The new physics 
introduced by Xi is mixing, with new dimensionless parameters, the species Schmidt 
numbers fio/fii, and the associated length scales for diffusion, the Batchelor scales. 
See Monin-Yaglom [23J for further information. 

Similarly to Theorem 1.1, for each fixed fi = (/jlq, • • • , ///), //^ > 0, j = 0, 1, • • • ,1, 
there exist a global periodic solution (u M ,xf, • • • , Xj)(^> x) in x G M 3 with period 
T P and a corresponding pressure function p^it, x) of the Cauchy problem fl6.1l) - (l6.2|) 
such that the equations in (16.1 1) hold in the sense of distributions, and the following 
properties hold: 

G L 2 (0, T; H 1 ) n C w ([0, T}; L 2 ) n C([0, T]; L Sl ), 
d t u» G L^Tjtf- 1 ) + (L S2 (0,T; ly" 1 '^) n L q (0, T; L r )), 

G L 2 ((0,T) x Tp) + L S2 (0,T;L^), 
Vp M G L 2 (0, T; if" 1 ) + L 9 (0, T; L r ), 
0<xf(t,x)<l, z = l,2,...,J, 
xf GL 2 (0,T;F 1 )nC([0,T];L S2 ), 

9 t xf e L 2 (0, T; H~ l ) + (L S2 (0, T; VF" 1 '^) n L 9 (0, T; L 2r )) , 
where 1 < s\ < 2, 1 < s 2 < oo, 1 < q < 2, and r = 2 (2g-i) ; and, in addition, 

d/-\u»\ 2 ) + V ■ (u"(i|uf + P n) + Mo|Vuf - ^oA(^) < f ■ u", 
^(xr)+V-(u^(xr))<^A/3(xf) 
in the sense of distributions, for any C 2 -function (3(x), 0"(x) > 0- 



KOLMOGOROV'S THEORY AND INVISCID LIMIT 15 

Theorem 6.1. The pathwise Kolmogorov hypothesis, Assumption (K41w), implies 
the strong compactness in L 2 fl L q of the velocity sequence u^(t, x) ; for some q > 2, 
when the viscosity \x tends to zero. That is, there exist a function u e L 2 (lL q ([0, T) x 
Tp) and a subsequence (still denoted) u M (t, x) such that 

u?(t, x) -> u(t, x) in L 2 n L q ([0, T) x T P ) as /x ->■ 0. 

Furthermore, 

0<xf(x)<l, z = l,...,/, 

which implies the weak-star convergence sub sequentially (still denoted) xf(^, x ) i n 
L°° to some functions Xj(t,x),0 < Xi(^? x ) < 1, as /x — >■ 0. Moreover, the limit 
function (u, xi, • • • , X/)(A X ) ^s a weafc solution of (j6.ip - fl6.2p and m addition, 

[ |u(t,x)| 2 rfx < / |u (x)| 2 dx + 2 / / u(s, x) • f (s, x) dxds. (6.3) 

JT P </T P JO </T P 

Denoting the weak-star limit in terms of Young measures, with v\ x iae Young mea- 
sure corresponding to the uniform bounded sequence xf for each i = 1, ■ - • ,/ ; we 
nave 

^=(^ x ,/3( Xj ))= W *-lim/3( X n 
for any C 2 -function fix- If, addition, (3"(x) > 0, then 



dtP{Xi) + V • G0(xi)u) < (6.4) 

m £/ie sense of distributions. 

Proof. The strong convergence of u M (t, x) can be obtained by following the same 
argument as for Theorem 5.1, which implies 

u**(t, x) -> u(t, x) in L 2 n L q ([0, T)xT P ),q> 2. 

Since < xf < 1, there exists a further subsequence (still denoted) xf an d an 
Z°°-function x« with < Xi < 1 such that 

Xi A X» as /x -»■ 0. 
Then the strong convergence of u M (t,x) implies that 

Xfu^(t,x)-x,u(^ x ) 

in the sense of distributions as /x — )■ 0. 
This implies 

+ V ■ (x»u) = 0, i = 1, • ■ • , /, 
in the sense of distributions. 

Furthermore, for any C 2 -function (3(x), P(Xi) * s uniformly bounded. By the 
representation theorem of weak limit by Young measures (cf. [26]; also see [I]), we 
conclude 

W) A Kx,/3(Xi)) as/x^O, 
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and, using the strong convergence of u M (t, x), we further have 

u^(t,x)/3( X f)-u(t,x)(^ x ,/3( Xi )) 

in the sense of distributions as fi — > 0. 
We now assume (3"(x) > 0. Since 

we take the limit — > above to conclude 

^Kx,/3te)) + V-(u(^ x ,/3( Xl )))<0 
in the sense of distributions. □ 

7. Implications to Numerical Computations 

We discuss what can be expected in terms of numerical convergence, on the basis of 
the convergence framework developed in this paper. Because the framework depends 
on compactness and subsequences, there is an implicit hint that uniqueness will be 
a problem. In fact, nonunique solutions for the Euler equations has been known for 
over a decade; also see Scheffer [25] and De Lellis-Szekelyhidi [HllTj. For numerical 
solutions, this means that LES solutions could depend in an essential manner on 
the SGS models or the equivalent in the form of modified numerical algorithms. 
Stated more directly, the turbulent Schmidt and Prandtl numbers of the models or 
the numerical Schmidt and Prantdl numbers of the algorithms could influence the 
selection of the numerical LES solution. Distinct but apparently converged solutions 
from distinct solutions for the identical problem are thus a possibility. Moreover, 
the LES solution regime is very sensitive to parameters (other than turbulent or 
numerical fluid transport) which introduce a regularizing length scale, and to the 
dimensionless ratios of these length scales. The literature for multiphase flow, tur- 
bulent mixing and turbulent chemistry has many dimensionless parameters. Those 
that are sensitive (experimentally, for example) in a particular problem can be ex- 
pected to be also sensitive in the definition of the SGS models and in their numerical 
manifestations as numerical analogues of the SGS models. Our optimistic hope is 
that a well designed LES algorithm should employ the SGS terms which eliminate 
further nonuniqueness due to numerical or algorithmic issues, and thus should yield 
a unique solution. 

Assuming uniqueness of solutions has been achieved, within an LES regime, what 
does convergence mean? Convergence of the velocity is strong (norm) in an L q 
space. This conventional notion of convergence requires no further discussion. For 
the passive scalars, the convergence is weak-star, meaning that convergence holds 
for any test function (observable) in the dual space (L°°)*. This dual space is a 
space of signed measures. The measures and the abstract space on which they are 
defined can be realized concretely in the form of the space-time dependent Young's 
measures. In language more familiar to physicists and applied mathematicians, 
the Young measures represent the space-time dependent probability distribution 
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function (pdf) of limiting values of the Xi- Although the limit also yields a classical 
weak solution, we expect that the pdf description will prove to be more useful. 
The reason for this claim is that nonlinear functions of the solution can pass to 
the limit in the pdf- Young measure formulation but not in the weak star limit 
formulation. The classical weak solution is the mean value of the Young measure 
solution, both regarded as space-time dependent distributions. The classical weak 
solution has far less information than the Young measure. Thus we state that the 
limiting concentrations Xi are pdfs. Clearly, this is a modification of standard ways 
of judging convergence, but now with a basis in mathematical theory. 

Is the Young measure limit a result of a weak mathematical technology, to be 
improved at some future date? Within the generality of Theorem 6.1, this seems 
unlikely. That theorem allows arbitrary and even infinite Schmidt numbers, meaning 
that sharp concentration discontinuities can persist for all time. Reasoning theo- 
retically, the essence of the Kolmogorov theory is that smaller and smaller vortices 
persist at all length scales, in the absence of viscosity. Thus the sharp interfaces 
must become ever more convoluted. Numerical evidence [TTl [18] supports this view, 
in which unregularized simulations display turbulent mixing concentration interfaces 
proportional to the available mesh degrees of freedom, uniformly under mesh refine- 
ment. More useful is to ask what can be expected in the case of finite, bounded 
Schmidt numbers. According to an analysis in |23j . the species diffusion has a loga- 
rithmetic singularity at high Schmidt number. But this argument omits time scales. 
A theoretical model [18], which does include time scales as well as SGS terms to 
modify the equations, suggests a rapid or instantaneous smoothing of species con- 
centration discontinuities, but at length scales that still may lead to Young measure 
solutions, i.e., convergence to a pdf. In view of the importance of the nature of this 
convergence, further studies, both numerical and theoretical, are called for. 

Finally, we note that our convergence framework stops at the level of passive scalar 
fields. Compressible Euler existence theories beyond this level are known only at 
the level of measure- valued solutions and even there at the level of convergent subse- 
quences, with no proof that the limit is actually a solution of the original equations. 
Because the concentrations are already described by this measure- valued framework, 
it is plausible that such a radical change in the convergence framework is required 
to move from passive scalars to true multiphase flow. Beyond the physical issues 
involved in the extension of Assumption (KJ^lw) to variable and species dependent 
viscosity and density, the mathematical issue is the proof of uniform estimates for 
partial differential equations whose coefficients are converging to Young measures 
(pdfs). 
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